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Abstract. We examine the role of anisotropic turbulence 
on the shear instabilities in a stratified flow. Such turbu- 
lence is expected to occur in the radiative interiors of stars, 
due to their differential rotation and their strong stratifi- 
cation, and the turbulent transport associated with it will 
be much stronger in the horizontal than in the vertical 
direction. It will thus weaken the restoring force which is 
caused by the gradient of mean molecular weight {/i). 

We find that the critical shear which is able to over- 
come the /i-gradient is substantially reduced by this 
anisotropic turbulence, and we derive an estimate for the 
resulting turbulent diffusivity in the vertical direction. 
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1. Introduction 

The recent observations performed by Herrero et al. (1992) 
reveal significant abundances anomalies in O and B stars. 
Even though these observations do not involve a high num- 
ber of stars, they display a striking correlation between the 
rotation velocity and the helium and nitrogen overabun- 
dances. 

Of all instabilities related with the (differential) ro- 
tation which have been examined so far, shear instabil- 
ities are the most efficient. They operate on a dynam- 
ical timescale, unlike the GSF instability (Goldreich & 
Schubert 1967 and Fricke 1968) or the ABCD instability 
(Knobloch & Spruit 1983), which both are of double dif- 
fusive nature and act on a thermal, hence much longer, 
timescale. In stellar radiation zones, the shear instability 
is hindered by the stable stratification, but nothing can 
prevent it if the flow is sheared in the horizontal direction. 
For this reason, as was argued by Zahn (1975), one should 
expect that the differential rotation in latitude induced by 
the meridian circulation maintains a turbulence which is 
strongly anisotropic, with a horizontal component of the 



turbulent diffusivity (Dh) much larger than the vertical 
one (D^). 

This turbulence tends to reduce the differential rota- 
tion which causes the instability, but it has little effect on 
the transport in the vertical direction. To obtain such a 
vertical transport, we have to invoke the instability gener- 
ated by the vertical shear, although it is impeded by the 
stratification. Fortunately, the stabilizing effect of the tem- 
perature stratification is weakened by radiative losses, as 
was first pointed out by Townsend (1958), and this mech- 
anism plays an important role in stellar interiors (Zahn 
1974). 

However in stars the stratification is due not only to 
the temperature gradient, but also to the variation of 
the mean molecular weight with depth, in particular out- 
side the regressing convective core of massive stars. These 
molecular weight gradients could well suppress the shear 
instability, if it were not for the anisotropic turbulence 
mentioned above. The horizontal diffusion associated with 
it will diminish the effect of the vertical molecular weight 
gradient, much like the thermal diffusivity reduces the ef- 
fect of temperature stratification. We shall examine in the 
present paper how this modifies the classical Richardson 
criterion and we shall give an estimate for the vertical 
diffusion coefficient. 

2. The effect of thermal diffusion on the shear in- 
stability 

The condition for the shear instability to occur may be 
obtained by stating that the work done against gravity to 
exchange two "bubbles" at levels z and z + Sz with veloci- 
ties U and U + SU must be smaller than the kinetic energy 
released by equalling their velocities (c/. Chandrasekhar 
1961); this yields the familiar Richardson criterion 



9 dAp/dz ^ 1 
PidU/dzf 4' 



(1) 
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where z is the vertical coordinate, g the gravity, p the 
density, and Ap is the density difference between the bub- 
ble and the exterior. It remains to evaluate that density 
gradient dAp/dz. 
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Quite generally, we may write 



d Alnp 



= 5 



dlnT dlnT' 



dz 



dln/i 
d^ 



Up 



[5 (V - V) + </.V^] , 



where the ' designates the conditions inside of the bubble, 
the gradients being expressed from now on with respect 
to the pressure, as is normally done in stellar structure 
theory: 



V = 



d InT 



and 



d In/i 



d InP 

We take the usual notations for 



(5 = - 



91np 
a InT 



and 



91np' 
9 In yLt / p y 



(2) 



(3) 



which depend on the equation of state (5 — X^^p = \ for a 
perfect gas). 

In the absence of thermal difFusion, V' is equal to the 
adiabatic gradient Vad and one retrieves the Richardson 
criterion relevant for a compressible fluid: 



< Ric 



where = + Nf, 



{dU/dz 

having split the Brunt- Vaisala frequency into 



tip 



and Nt 



(4) 



(5) 



with Ric ~ 1/4 being the critical Richardson number. 

The effect of thermal damping has been examined by 
Dudis (1974), and Zahn (1974) has given an estimate for 
the turbulent diffusivity in the limit of low Peclet num- 
ber {vljK <c; 1, with K being the thermal diffusivity). 
Recently, Maeder (1995) derived a new expression for this 
coefficient, which has the advantage of describing the in- 
termediate regimes. It is inspired by the mixing length 
treatment of non-adiabatic convection (see Cox & Giuli 
1968), where one defines the ratio F between the thermal 
energy transported by an eddy of volume V 



pVCpAT = pVC,,T—{V - V) 
tip 



(6) 



and the energy which diffuses through its surface of area 
A during its lifetime ijv 



\-w) 



tip 



(7) 



(AT is the maximum temperature difference, AS* is the 
entropy excess, and x = pCpK ~ WaT^/Snp the ra- 
diative conductivity, with the usual notations). Assuming 



that the eddy may be described as a sphere of diameter i 
{i.e. V/A = one has 



r = 



6K 



V - V 

Vad - V 



(8) 



This yields the modified version of the Richardson cri- 
terion which was established by Maeder (1995): 



9- 



dAlnp 

dz 



r + 1 



dz J 



(9) 



According to this criterion, the vertical shear can be un- 
stable only provided 



dU 
d7 



> 



1 



(10) 



Then eq.(||) may be solved for the parameter V = vH/QK 
which characterizes the largest eddies that satisfy the 
Richardson criterion, and these produce a turbulent diffu- 
sivity = 2KT. 

Recently Meynet & Maeder (1996) calculated the evo- 
lution of rotating massive stars with this prescription for 
the turbulent difFusion, assuming local conservation of 
angular momentum. They found that, even though the 
model displays a strong differential rotation at the end of 
the main sequence, no significant mixing occurs because 
the molecular weight gradient inhibits the shear instabil- 
ity. Since according to the observations such mixing does 
take place in fast rotators, they conclude that this pre- 
scription underestimates the shear mixing in stars, and 
there must be some way to overcome the /i-barrier. 

3. The effect of the horizontal diffusion 

As was already mentioned in the Introduction, it seems 
plausible that stars exhibit an anisotropic turbulence due 
to the horizontal shear which is induced by the meridian 
circulation (see also Zahn 1992). This turbulence gener- 
ates a horizontal transport which is much larger than the 
vertical one, and which tends to smooth out the horizon- 
tal fluctuations of molecular weight. Therefore p will vary 
along the trajectory of the bubble considered above (see 
fig-1). 

The restoring force per unit mass is then given by 

with the same notations as before. 

To estimate (V^ — V^^) we proceed as above for the 
thermal difFusion of heat: it suffices to replace AT by A/x 
and K by the turbulent diffusivity D^. The result is 



V' — V — - V 



where now — v£/6Dh 



(12) 
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Radiative losses 



Horizontal diffusion 



It is an easy matter to solve this second order equation 
15] ) , and to derive from it the value of vertical diffusivity 
Dy ~ vi/3. Solutions exist only if 



dU 
d7 



> 



QRic 



Ar2 



K + Dh Dh 



(16) 



which replaces the much more stringent condition (pO|). 

Note that the vertical turbulence transports heat and 
thus modifies the temperature gradient V. This was taken 
into account by Maeder & Meynet (1996) when they evalu- 
ated the mixing by shears. However, in most circumstances 
relevant to radiative interiors, we are in the presence of 
mild turbulence, so that x <C K, and we may ignore this 
effect. Furthermore, we stated as a working hypothesis 
the condition Dy ~ 2a; ^ Dh- The turbulent diffusivity is 
then given by the simple expression 



Dy 



2RicidU/dz) 



Nll{K + Dh) + miDh 



(17) 



Fig. 1. As a turbulent eddy moves vertically, it is partially 
mixed with the surrounding medium, and therefore its compo- 
sition varies. 



We thereby obtain still another version of the Richard- 
son criterion: 

r 



r + 1 



T 



r, 



dC/ 
d^ 



(13) 



which describes how the stabilizing effect of the /i-gradient 
is diminished by the horizontal turbulence. Note that the 
influence of the temperature gradient will also be modified 
by this turbulence since the advective transport of heat 
will be added to the radiative diffusion: K ^ [K + Dh)- 
We see that any vertical shear may become unstable, 
for small enough vi- But the Reynolds number character- 
izing these eddies must be larger than some critical value, 
otherwise their motion will be damped out through vis- 
cous friction. Hence they must obey the condition 

— < Rec, (14) 

where v is the molecular viscosity, and Rec ~ 10 is the 
critical Reynolds number. 

Among all turbulent eddies satisfying (13) and (p^, 
those contributing most to the vertical transport arc those 
for which vi is the largest. Their size will be given by the 
equality in eq. (|l3|), namely 



One final word about the value of the ratio V/A be- 
tween volume and area of the turbulent eddies. Its choice is 
rather arbitrary. Following Maeder (1995), we have taken 
V/A — £/6, whereas Bohm-Vitense (1958) proposed the 
value V/A = (2/9)^ when she developed her famous mix- 
ing length formalism. We could refine this numerical coef- 
ficient, for instance by taking our inspiration from Dudis' 
(1974) work, where T = {3/8)v£/K- However, we prefer 
to hide this uncertainty in the choice of Ric and in the 
prescription we use to evaluate the horizontal diffusion 
coefficient Dh , which also involves an unknown parameter 
of order unity (see Zahn 1992). 

4. Conclusions 



We have shown how the Richardson criterion, which gov- 
erns the shear instabilities in a stratified fluid, is modified 
in the presence of a strongly anisotropic turbulence. This 
turbulent erosion is similar to that which intcrfers with 
the advective transport of a large-scale flow, and which 
has been described by Chaboyer and Zahn (1992). An- 
other example of this phenomenon is given in a numeri- 
cal simulation performed by Vincent et al. (1996). Here 
it diminishes the stabilizing action of a molecular weight 
gradient, to a point which makes it plausible that rapidly 
rotating stars be partially mixed. The prescription (17) 
we have derived above for the vertical diffusivity is being 
implemented in our stellar evolution codes, and the results 
will be reported in a forthcoming paper. 



-m 



x + K + Dh' ^ x + D'"'' 



where x = u£/6 > vRcc/Q- 



Nf, = Ric 



dU 



(15) 
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